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Abstract 

We develop a method for systematically constructing Lagrangian functions for dissipative mechanical, 
electrical and, mechatronic systems. We derive the equations of motion for some typical mechatronic 
systems using deterministic principles that are strictly variational. We do not use any ad hoc features 
that are added on after the analysis has been completed, such as the Rayleigh dissipation function. 



We generalise the concept of potential, and define generalised potentials for dissipative lumped system 
elements. Our innovation offers a unified approach to the analysis of mechatronic systems where there 
are energy and power terms in both the mechanical and electrical parts of the system. Using our novel 
technique, we can take advantage of the analytic approach from mechanics, and we can apply these pow- 
erful analytical methods to electrical and to mechatronic systems. We can analyse systems that include 
non-conservative forces. Our methodology is deterministic and does does require any special intuition, 
" 53 ■ an d is thus suitable for automation via a computer-based algebra package. 

a. 

Introduction and motivation 

It is a widely believed that the Lagrangian approach to dynamical systems cannot be applied to dissipative 
systems that include non-conservative forces. For example, Fcynman [1] writes that "The principle of least 
action only works for conservative systems — where all the forces can be gotten from a potential function. " 
Lanczos [2], writes "Forces of a frictional nature, which have no work function, are outside the realm 
of variational principles, while the Newtonian scheme has no difficulty in including them. Such forces 
originate from inter-molecular phenomena, which are neglected in the macroscopic description of motion. 
If the macroscopic parameters of a mechanical system are completed by the addition of microscopic pa- 
rameters, forces not derivable from a work function would in all probability not occur. " Lanczos [2], also 
writes "Frictional forces (viscosity) which originate from a transfer of macroscopic into microscopic mo- 
tions demand an increase in the number of degrees of freedom and the application of statistical principles. 
They are automatically beyond the macroscopic variational treatment. " It is clear that dissipative forces 
present a problem to traditional Lagrangian analysis, which means that the Newtonian approach has 
historically had an advantage, particularly where dissipative forces are significant. 

There are a number of formalisms for applying a Newtonian (force-based) approach to mixed mecha- 
tronic systems. The bond-graph approach is based on the systematic use of effort and flow variables. The 
work of Karnopp et al. [3] is important in this regard. We will employ some aspects of Karnopp's work, 
including the homomorphic mappings of variables between different systems. There are clear analogies 
between mechanical and electrical oscillators, and we make use of these. 



The Newtonian approach has been dominant in practical discipline areas, such as mechanical engi- 
neering. In contrast, the Lagrangian approach, which is very elegant, has tended to dominate advanced 
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physics texts. For example, the Hamiltonian approach dominates the subject of quantum mechanics. 
Penrose [3] , refers to this paradigm as the "magical Lagrangian formalism. " He goes on to write that 
"The existence of such a mathematically elegant unifying picture appears to be telling us something deep 
about our physical universe. " 

There are a number of more prosaic factors in favour of the Lagrangian approach, which include: 

• In the Lagrangian formulation, forces of constraint do no work and need not be considered in the 
analysis. It is often not necessary to calculate internal stresses or forces of reaction. 

• Post [5] points out that it is easy to state the underlying physical laws in arbitrary, curvilinear 
coordinates. It is possible to use generalised coordinates that directly reflect the nature of the 
physical system. 

• Noether's theorem tells us that, if the Lagrangian function possesses a continuous smooth symmetry, 
then there will be a conservation law associated with that symmetry [4]. For conservative systems, 
this leads to the laws of conservation of momentum and conservation of energy. These conservation 
laws essentially give us one integration of the laws of motion for free. For example we can calculate 
the final momentum, and the final energy of a system without the need to explicitly integrate the 
laws of motion. 

• Lagrangian modelling of machines, automatically takes care of energy transfer between different 
components of a whole system. This prevents incomplete models, which give rise to errors and 
paradoxes, such as the problem of the Penfield motor [6]. We believe that Lagrangian modelling 
is a natural choice, where energy is exchanged between different types of storage elements, in such 
systems as: a moving wire in a magnetic field, the D'Arsonval moving-coil meter, or for mechatronic 
systems more generally. 

• In the Hamiltonian formulation, only first derivatives are required, not second derivatives. 

• Many quantum systems, such as the hydrogen atom, only have a few degrees of freedom and a 
complete description of all the microscopic parameters is possible. This means that frictional forces 
may not even need to be considered. 

Perhaps the strongest theoretical motivation for the Lagrangian approach is that it explicitly rep- 
resents the symmetries of the underlying physical laws. Melia [7] writes: "As we shall see, the sole 
motivation for using action principles is to improve our understanding of the underlying physics, with a 
goal of extracting additional physical laws that might not otherwise be apparent. " 

Prior to the work of Riewe 0|9], there was no satisfactory method for completely including non- 
conservative forces into a variational framework. Riewe writes that "It is a strange paradox that the most 
advanced methods of classical mechanics deal only with conservative systems, while almost all classical 
processes observed in the physical world are non-conservative. " We regard the approach used by Riewe 
as the most satisfactory method for including non-conservative forces into a variational framework. In 
this paper we apply his approach, for mechanical systems, to the new areas of electrical and mechatronic 
systems. We also provide recipes for constructing Lagrangian functions and show, by example, how these 
techniques can be employed more generally. We believe that the Lagrangian approach naturally models 
energy exchange within complex machinery, where energy can be stored and transferred between many 
different forms, including: energy of inertia, elastic energy, frictional loss, energy of the magnetic field, 
energy of the electric field and resistive loss. Our approach can be used to confer the advantages of the 
variational method of analysis to a wide range of mechatronic systems, including systems that suffer from 
dissipative loss. 
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A short summary of the variational approach 

We can denote a Lagrangian function for a system as C, then we can specify the total action of the system 
as ^ 

1= [ 2 Cdt, (1) 

where T\ and T2 represent the boundaries of the closed time interval over which we wish to conduct 
our analysis. Equation [T] is referred to an action integral. It is a functional that maps functions, C, 
onto numbers, T. The Euler-Lagrange equation specifies necessary conditions for the first variation of 
the action integral to vanish, 8 [I] = 0. Suppose that the Lagrangian function includes references to a 
generalised coordinate, x(t), and to its first derivative x so C = C (x, x), then the action is extremal when 
we choose x(t) in such a way that the Euler-Lagrange equation is satisfied: 

d fdC\ dC n 

dt{dx)-0-x-°- ® 

This is the same as saying that all first order variation of the action is zero, 8 [1] = 0. The Eulcr Lagrange 
equation is an ordinary differential equation that describes the dynamics of the system, in terms of the 
specified generalised coordinates, such as x{t). 

For mechanical systems the Lagrangian is written in terms of energy functions, which are summed 
together with appropriate sign conventions. They typical symbols are kinetic energy of inertia, T (£), and 
potential (elastic or gravitational) energy, V (x). For these systems the Lagrangian function can be written 
as: £ = T—V. As we shall see, a classical example is a mass on a spring, where C = T—V = \mx 2 — ^kx 2 . 

We will use the notation of Gel'fand [TU], who denotes a general fcth order derivative as: x^- k \ This 
is more versatile than the more traditional "dot" notation, of Newton. It is common for Lagrangian 
analysis to only consider integral derivatives, of low orders, of the generalised coordinates. For example, 
we might consider x = x^°\ x = x^ and possibly 'x = x^ 2 \ Gel'fand writes the generalised integer-order 
Euler-Lagrange equation a form that includes higher derivatives, and is equivalent to: 

e<-')'4(^t)= ' « 3 > 

fc=0 v 7 

where where k £ Z, and where it is understood that J^f k) = 0, for values of k where the Lagrangian 
has no dependence on the k'th derivative of the coordinate. The proof of Equation [3j can be obtained 
by repeatedly integrating by parts, and applying the du Bois-Rcymond lemma. Proofs can be found in 
Gel'fand [TU] and Smith [UJ. 

Since the seminal work of Riewe [H1H], a number of other authors have used his approach. These 
include Agrawal [12], Rabei [13], Frederico [14], Musielak [15], Elnabulsi [16] and Almeida [TT1[T8], but 
none of these papers have addressed the problem of dissipative forces in electrical or mechatronic systems, 
which is our main concern here. 



Fractional Calculus 

The indices of differentiation in The Euler Lagrange Equation [3] can be fractional, which leads to the 
formulation: 

\ ^ d k f dC \ 



^ ■5?lfaWj= ' (4) 
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where k £ Q, and where it is understood that d ^f k) = 0, for values of k where the Lagrangian has no 
dependence on the /c'th derivative of the coordinate. The proof of this proposition depends on a fractional 
version of integration by parts, and is found in Riewe [5]. 

The history of the fractional calculus can be traced back, at least, as far as a paper by Liouville in 
1832, where he gives a definition of a fractional derivative: If a function, f(t) can be written as a sum of 
exponential functions, or possibly complex-exponential functions, as 

/(t) = ^c m e a - 4 , (5) 

m 

for a suitable choice of values, c m and a m , then we can write the fractional derivative as: 

fW d k f{t) _d k ^ 



d(t + oo) k d(t + oo) 



k / < 



E 1 



l (o m )*e°» t . (6) 



Fractional derivatives are not local if k is not an integer, which means that their value depends on a 
region around the point of evaluation. The choice of region is important. There is also the question as 
to whether the derivative is causal, or not. Riemann defined the fractional derivative as: 

,f vf ® = =^ ■ [ (t- t'y^f^dt' (7) 

J d{t-c)~v T(v) J Tl 

provided that $t(v) > 0. We argue that the derivative is causal if it is a function of the state of the system 
over the time interval that includes all time up to the present moment, — oo < t' < t. It can be shown that 
Riemann's definition is consistent with Liouville's 1832 definition, provided T\ — > — oo, which means that 
Liouville's definition is causal. In this paper we use Liouville's definition of fractional derivatives. This 
means that we have to express any functions of interest as sums of exponential functions as indicated in 
Equation [5] The responses of systems can often be written in this form using a finite number of terms, so 
Liouville's definition is very convenient. In the more general case, we can use the theory of the Laplace 
transform to perform the necessary decomposition. In particular, we can write: 

^M=L^{s^L{f}(s)}. (8) 

Once again, care has to be taken to integrate the Laplace transform over the correct time interval. This 
can be achieved by using step-functions as windowing functions, and by making use of shifting theorems. 
The theory of the fractional calculus has been well documented, and summaries can be found in Old- 
ham et al. []j|] and briefly in Wiesstein [2"0] . 



A mechanical harmonic oscillator 

We consider a common problem from classical mechanics, of a mass on a spring. This problem is widely 
used to define notation, and can be found in: Lamb [21], Goldstein [22], McCuskey [23], Resnick 
& Halliday [24], Whylie [25], Fowles [26], Feynman [TJ, Rabenstein [27] and Lomen [28], and many others. 

The mechanical harmonic oscillator consists of a mass, spring and massive support (or foundation). 
The complete system is shown in the schematic diagram in Figure [T] A mass, to, is attached to a springj, 

lr There is some difficulty with the schematic notation for the spring, k, since the traditional schematic symbol for a 
spring resembles the traditional schematic symbol for a resistor. This creates problems if we need to represent both of these 
different objects in a single drawing. We have followed examples from Giesecke et al. 1291 . In particular our symbol for the 
spring has a different aspect ratio to the symbol for the resistor, and the terminations at the ends are different. 
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Figure 1. A mass on a spring: We consider the simple introductory problem of a mass, m, on a 
spring, with Hookc's-law constant of k. The kinetic energy stored by the inertia of the mass is denoted 
by T = \mx 2 . The elastic potential energy stored in the spring is denoted by V = \kx 2 . The 
independent coordinate is denoted by the position, x. The Lagrangian function is traditionally written 
as C = T — V, which can be written explicitly as C = ^mx 2 — \kx 2 . 



k, which is attached to a massive support. The position of the spring is measured relative to a datum 
position, which is in a fixed position relative to the massive support. Without any loss of generality we 
can choose the location of the no-load position of the mass, which gives a simple rule for the stored energy 
in the spring, V = ^kx 2 . 

The classical problem of a mechanical oscillator is shown in Figure Q] Together the mass and spring 
form a mechanical harmonic oscillator. Williams [30] specifies a Lagrangian for this physical system in 
terms of the single spatial coordinate, x, and writes: 

£ = T-V = -mx 2 --kx 2 7 (9) 
2 2 w 

where T = \m (i) 2 , which is the kinetic energy in the inertia of the mass and V — ^k (x) 2 is the strain 
energy stored in the spring. The Lagrangian function, in Equation [H] is in a form where we can directly 
apply the procedure of Euler and Lagrange Equation to obtain: 

mx + kx = mx {1) + kx {0) = 0, (10) 

which is the standard Ordinary Differential Equation (ODE) for this system. Equation [TT] can be solved 
using a number of techniques, including the method of the Laplace transform, to obtain: 

x(t) = A cos (cot) + B sin (cot) , (11) 

where A and B should be chosen in order to satisfy the initial conditions and, uj = yjk/m is the un- 
damped natural angular frequency of oscillation, in radians per second. In the completely un-damped 
case, w is also the resonant angular frequency. We note that the solutions for x(t) can be written as 
sums of exponential functions, since cos (#) = fe +J + e"-? 6 ') /2 and sin (6) = (e +: > — e~ 3 1 ) /(2j). This 
means that Liouvillc's definition in Equation [5] could be directly applied if we wanted to determine any 
fractional derivatives. 



A homomorphic mapping 

The example shown in Figure[T]is simple and well known, and lies completely within a mechanical problem 
domain. It is not immediately obvious how to extend this type of work to an electrical domain. We need 
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a homomorphic mapping of variables that can relate different variables in different physical domains The 
mapping needs to relate the names of variables, as well as the set of permissible functions and operators, 
that work on those variables. We use the mapping described in Karnopp et al. [3], which is summarised 
in Table [U 

Table 1. Table [TJ A homomorphic mapping due to Karnopp et al. 





1YT ppb a n i pa 1 

. 1 1 V . 1 11(1111 1(11 


"Fjlpftripal 


coordinate 


displacement, x 


charge, q 


flow variable 


velocity, v = ^jf 


current = i = 4| 


energy 


energy U 


energy U 


effort variable 


force, F = % 


voltagc= v = ^ 


energy increment 


dU = F ■ dx 


dU = v ■ dq 


power = dlA / dt 


<p = F-^ = F-v 


V = v ■ — v ■ i 


inertial element 


mass, m 


inductance, L 


generalised momentum 


momentum, p = m ■ 


magnetic flux, $ = L ■ $ 


Newton's second law 


F = m ■ 4t# = m ■ a 


„ _ r, . £q — T , di 


elastic element 


Hooke's law constant, k 


inverse capacitance, 1/C 


Hookc's second law 


F = —k ■ x 


v = +7J ' 1 


dissipative element 


damping, c 


resistance, R 


frictional force 


F = -c-& = -C'V 

rt dt 


V = +R- % = +R-i 


Joule's law 


V = v l c 


V = i A R 


energy of inertia 




U= if = ±L-% 1 


clastic energy 







A homomorphic mapping: The names and purposes of the most important mechatronic dynamical 
concepts. 



An electrical harmonic system 

If we place a capacitor, C in parallel (and series) with an inductor, L, as shown in Figure [U then the 
resulting system will form an electromagnetic harmonic oscillator. 

If we temporarily ignore the presence of resistance, then we obtain the circuit in Figure [5] is the 
exact analogue of the mechanical system in Figure [T] In order to emphasise the homomorphic mapping 
between the mechanical and electrical domains, we map the Lagrangian function in Equation [3J using 
the homomorphic mapping in Table [T] to obtain 

where L is the inductance, C is the capacitance and q is the charge that is transferred through the circuit. 
We note that Equation [12] is a correct Lagrangian function for the circuit in Figure [2j 

In Equation lI21 we use the variable q as a coordinate, in accordance with homomorphic mapping due 
to Karnopp. It is the more usual practice in electrical engineering to use the voltage across a capacitor, v, 
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Figure 2. An LC electromagnetic harmonic circuit: We consider a capacitor, C, in parallel with 
an inductor, L. We consider the idealised case where there is no dissipative loss, or resistance R. The 
Lagrangian function can be written as C = T — V, where the magnetic energy stored by the field of the 
inductor, denoted by T = L q 2 /2, and the electrical potential energy stored in the capacitor is denoted 
by V = q 2 / (2C), and q = Cv is the coordinate, which we interpret as the electrical charge that is 
transferred through the circuit. 



as though it were a generalised coordinate. Fortunately, it is possible to subject the Lagrangian function 
in Equation [12] to a Legendre transformation, of q = Cv, to obtain a new Lagrangian function that is 
consistent with the previous Lagrangian function (in terms of energy exchange) , but uses the conventional 
coordinate of v (rather than q). This new Lagrangian function is self-contained in the sense that the 
energy terms for inductor only include references to parameters that pertain to the inductor, and the 
energy terms for capacitor only include references to parameters that pertain to the capacitor. There 
are no cross-terms. If we impose this last condition then the Legendre transformation is unique and we 
obtain a new Lagrangian function: 



C=T-V= 




where the independent generalised coordinate is now v. The function is the zeroth derivative of v, 
which is identical with v. We can write i/ ' = v. The function denotes the derivative of v to the 

order of -1, which is equivalent to the integral of v. We can write i/" 1 ) (<) = J T v (r) dr. 



Lagrangian terms for some common lumped mechatronic elements 

We can see from the last example that electrical and mechanical systems can be mapped onto one an- 
other but some care has to be taken with regard to what we regard as a coordinate. The canonical choice 
for a massive particle is to regard the spatial position as the coordinate and to regard the generalised 
momentum as the other variable of interest. These choices are not arbitrary. The coordinate, x, must 
be an exact differential, for example: § xdx = for all possible closed paths. In quantum mechanics, 
position, x, and momentum, p, arc conjugate variables. The relationship between x and p is a physical 
phenomenon, not just an arbitrary choice. Finally we know from classical mechanics that x and p play 
a role in Liouville's theorem. See Reif [31] and Wannier [32], for example. Liouville's theorem would not 
apply in the same way if we were to describe particle motion in terms of force and velocity, rather than 
position and momentum. If we want to apply Liouville's theorem to a complicated electrical system with 
many degrees of freedom then we really need to use magnetic flux, and electric flux, q, to describe 



8 



each element of the system. This has been carried out in some specialist areas, such as Allison |33| . 
but it is not common and is not likely to become universal in the electrical engineering literature in the 
foreseeable future. 

The conventional choice of electrical variables, voltage, v, and current, i are of a fundamentally differ- 
ent type to the conventional choice of mechanical variables, x and p. These incompatible conventions arc 
not likely to change. The best solution seems to be that we should re-write the Lagrangian terms, using 
the conventional electrical variables, but to do so in such a way that the energy values are preserved, and 
all sign and phase relationships are consistent. These conversions are not physical laws. They have the 
same status as the conversion from degrees Fahrenheit to Kelvin, for example. The Lagrangian terms, 
using common electrical engineering variables, are shown in Figures [2] and [3J If we wish to include one of 
these devices in a system, and we wish to carry out Lagrangian analysis then it is only necessary to look 
up the term in the table and to include it in the appropriate way, in the Lagrangian function. No other 
special modelling or cognition is required. The process is direct enough to be able to be carried out by 
computer. 



Table 2. Table of Lagrangian terms, in terms of current 



parameter 


phasor 


Lagrangian C 


order, k 


Euler-Lagrange 


L 


I=i-j/(uL))V 









R 


I=(1/R)V 


+2-R(it-V*)f 


-1/2 




C 


I = (+juC) V 




-1 





Lagrangian terms, with current: We list the common electrical lumped parameters, and compare 
the admittance with the corresponding term from the Lagrangian function. We also list the order of 
differentiation, k, and the corresponding term from the Euler-Lagrange equation. The phase direction 
of the Lagrangian term leads the phase direction of the admittance by 90°. This is equivalent to 
multiplying the Lagrangian term by +j. We can multiply the Lagrangian term by any constant that we 
like, as long as we do this consistently. If we were to remove the factor of +j then the Lagrangian terms 
and the admittances would have consistent phases, but all the Lagrangians would have new phases and 
some of these would not be consistent with existing practice in mechanics. In this paper, we rigorously 
adopt the convention that is used in mechanics, which means that we do not use the sign convention 
that is common in electrical engineering. 

In the interests of consistency, we include a table for the more familiar Lagrangian terms, in Table 01 

Finally, we note that some of the Lagrangian terms are imaginary and that the resulting Lagrangian 
functions will, in general, be complex. We consider arithmetic operations to operate within the field of 
complex numbers, C. The traditional case, of real Lagrangian functions, is a special case of our more 
general formulation. Our formulation is consistent with the earlier work of Illert [34] , who applied the 
concept of complex Lagrangian functions to the classical seashell problem. 



A damped mechanical harmonic system 

We consider the damped mechanical oscillator, as shown in Figure[31 This problem is solved by Riewe [3]. 
In our case, we only need to take the Lagrangian terms from Table 0] and add them together to form the 
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Table 3. Table of Lagrangian terms, in terms of voltage, v; Table [3] 



parameter 


Phasor 


Lagrangian C 


order, k 


Euler-Lagrange 


L 


V=(jujL)I 




-1 


-l W (-2) 


R 


V = RI 




-1/2 


—M-i) 


C 


V = -juCI 










Lagrangian terms, with voltage: Wc list common lumped electrical parameters, and compare the 
impedance with the corresponding term from the Lagrangian function. We also list the order of 
differentiation, k, and the corresponding term from the Euler-Lagrange equation. The phase direction 
of the Lagrangian term lags the phase direction of the impedance by 90° . This is equivalent to 
multiplying the Lagrangian term by —j. We can multiply the Lagrangian term by any constant that we 
like, as long as we do this consistently. If we were to remove the factor of — j then the Lagrangian terms 
and the admittances would have consistent phases, but all the Lagrangians would have new phases and 
some of these would not be consistent with existing practice in mechanics. In this paper, we rigorously 
adopt the convention that is used in mechanics, which means that we do not use the sign convention 
that is common in electrical engineering. 



Table 4. Table of mechanical Lagrangian terms, Table 0) 



parameter 


Lagrangian 3D 


Lagrangian ID 


order, k 


Euler-Lagrange 


m 


1 1 |2 




1 


—mx^ 


c 


not applicable 




1/2 


-cx« 


k 


-ifc|x| 2 


-iJb(*o») J 








Lagrangian terms for mechanical parameters: We list the common lumped mechanical 
parameters. From left to right, we list the common symbol for the parameter, the Lagrangian for the 
3D vector case (in terms of position or momentum) , the Lagrangian for the ID case (in terms of the 
position only), the order of differentiation employed and the resulting them in the Euler-Lagrange 
Equation. We use the sign convention of Gel'fand and Fomin [10] for the terms of the Euler-Lagrange 
Equation. Some authors introduce an additional minus sign, in order to make all of these terms positive. 

Lagrangian function for the system. This is shown as follows, 

C = \m{z)* + ic(xW*>y -\k{xf 

= I^^W^ + ^^d^^-Ifc^W) 2 . (14) 

The resulting Euler-Lagrange equation can be assembled from the Euler-Lagrange terms in Table 21 or 
calculated directly, using Equation SJ The result is given by 

- (m.T (2) + cx (1) + fcr (0) ) = 0, (15) 

which is the same result that we would obtain if we used a free body diagram and Newton's laws of motion. 
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Figure 3. A damped mechanical harmonic oscillator: By introducing non-conservative, or 
dissipative, elements into a system we need to generalise our concept of potential. The Lagrangian for 
this system can be written as C — \vn (a^ 1 )) — |fc (a; 1 - - 1 ) + |c (a^ 1 / 2 ') . The additional term differs 
from the terms for the un-damped system in two key ways: the term is complex has an imaginary 
phase, of +j, and there is a fractional derivative of the coordinate, a^ 1 / 2 ). 



The use of constraints 

The use of the calculus of variations to evaluate extremal functions, subject to constraints is described 
in a number of references, including Lanczos j2j. It is possible to regard perfect sources, of voltage or 
charge (or force or velocity), as constraints. This can simplify the working of some problems. We include 
an example here. 

We consider the case of a purely resistive system. Jaynes [35] traces this problem back as far as 
Kirchhoff [36| , and points out that the condition that no electric charge should accumulate at any point 
in a resistive material requires that V • (aVv(x)) = 0, which is just the Eulcr-Lagrange equation stating 
that the production of Joule heat in a domain T>, J v a (S7v) dT> is stationary with respect to variations 
5v (x) that vanish at the boundary of T>. 

We also note that Kirchhoff's voltage law is a tautology based on the definitions of energy and volt- 
age. In a quasi-static situation where radiation is not significant voltage is just energy per unit charge, 
v = AU/Sq. It is tempting to regard Kirchhoff's current law as a statement of the conservation of charge, 
but this is misleading. Even if we grant the continuity of charge, V J + dp/dt ~ then it would still be 
possible to have accumulation of charge. The equivalent principles of "no accumulation of charge" and 
"minimum production of Joule heat in a domain" are ultimately statistical in nature and are related to 
the second law of thermodynamics. This is discussed in Allison [53] , 

In this paper, we consider the special case where the parameters are lumped into two resistors in 
series, as shown in Figured] To simplify the notation, we define a gradient operator, V, over variations 
with respect to the variational operators, S Vl and S V2 , rather than partial derivatives, d/dv\ and d/dv2- 
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Wc can define 



V£ = 



5 Vl C 
S V2 C 



(16) 



This allows us to use the notation of the Lagrange multiplier to write down a necessary condition for a 
constrained optimum. If the independent variables, v\ and vi are constrained by a function of constraint 



vi 

V2 



0, 



then we can only obtain constrained stationary values of C when 

V£ - AW = 0. 
We apply this principle to the problem in the next section. 



(17) 



(18) 



A problem with two resistors 

In Figure IHa), the voltage source, vs places a constraint on the voltages across the two resistors, v\ and 
V2- We can apply Kirchhoff's voltage law to the single mesh in this circuit to obtain <!> = v\ +1)2 — vs = 0, 
where ^ is a function of constraint. Gradients of this constraint function are needed in order to determine 
the constrained stationary functions of the system. Kirchhoff noted that the voltages in a resistive circuit, 
such as Vi and vi , would arrange themselves in such a way as to minimise the dissipated heat energy, given 
by Joule's law. Jaynes points out that this is equivalent to defining a "Kirchhoff" Lagrangian function. 
For us, this takes the form Lk = power = v\ 2 / R\ + V2 2 / R.2- In order to ensure compatibility with other 
existing Lagrangian functions, we apply a Legendre transformation to obtain a new Lagrangian function 




The resulting Lagrangian function, in Equation 1191 can also be assembled from the Euler-Lagrange terms 
in Table [3J It is consistent. In Figure SJb) wc indicate that this situation is very common since it occurs 
whenever a linear source is connected to a resistive load. 



We can use the principle of the Lagrange multiplier to obtain equations for the stationary functions, 



subject to constraints. We begin by calculating the individual variations: S Vl C = —(1/R\)v 
, which leads to the following form for the gradient, V£, as 



S V2 C 



(l/R2)vi~ 1) 



V£ = 



5 Vl C 






5 V2 C _ 




L „„(-i) 



The function of constraint is 







' 1 " 


. s V2 y _ 




1 



Vi + t>2 — vs — 0, and wc obtain the gradient of this as 

r ; 

We can apply the principle of the Lagrange-multiplicr to obtain 

A 



Hi V1 
R 2 2 



0, 



(-1) 

1 



and 



(20) 



(21) 



(22) 
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Figure 4. Two resistors in series: The voltage source, vs places a constraint on the voltages across 
the two resistors, v\ and v 2 . Kirchhoff's voltage law implies that * = V\ + v 2 — v$ = 0. We can regard 
the function 'I' as a function of constraint. Kirchhoff used a function that is equivalent to the dissipated 
power as a Lagrangian function. In modern notation we can write Ck = P = /Ri + («2) jRi- In 
order to ensure compatibility with other existing Lagrangian functions, we need to apply a 
transformation to obtain a new Lagrangian function of C — t^t- (tV" 1 / 2 ') + (^2^ -1 ^ 2 ') ■ This 
Lagrangian function has been multiplied by a scalar of —j/2 and the order of differentiation has been 
reduced from to i>( -1 ' 2 ), which is equivalent to a half-integration of the Lagrangian function, or a 
full integration of the resulting Euler Lagrange equation. 



for some constant complex number, A G C. If we differentiate once, with respect to time, and multiply 
by —1, we obtain: 

* = £ = (23) 
i?l R 2 

where i € C is a complex number. Since i>i, v 2 and i?2 are all real we an infer that i is real. If we 
consider Ohm's law then i is a common current that is shared by both resistors. Equation [231 can also be 
obtained by using Kirchhoff's current law and by applying Ohm's law twice. Our main aims in presenting 
this last example are: 

• to illustrate the use of constraints, with possible time and rate dependence, 

• to demonstrate the utility of our extended Lagrange-multiplier notation, 

• to provide a historical reference to the important work of Kirchhoff and Jaynes, and 

• to resolve an apparent contradiction between the Lagrangian analysis of purely reactive systems 
which only have energy storage, and the Lagrangian analysis of purely resistive systems, which only 
have power dissipation. The reactive systems have Lagrangian terms that only depend on energy 
terms. The resistive systems have Lagrangian terms that only depend on power terms. Superficially, 
this appears to be a contradiction. 

In order to examine the proper relationship between "energy" Lagrangian terms and "power" La- 
grangian terms further, we next consider a mixed example, where the system contains both resistive and 
reactive elements, tied together with a constraint. 



A damped electrical harmonic system 

We can use the terms in Table [2] to write the Lagrangian function for the circuit in Figure [5] as: 

C^ + l«(,^)\ 1 -L(^>y-±(, c ^)\ (24, 
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We can use Kirchhoff 's current law to impose a constraint function of VP = in + zl + ic — is = 0- We can 
then use the techniques from the last section to obtain the solution to the problem of the dynamics of 
this circuit. We then find a functions, v(t), that give stationary values for the of the action, X, subject to 
the constraint, ty. The principle of the Lagrange multiplier allows us to replace the optimising principle 
with a new constraint, 

V„ [J] - A ■ V„ [*] = 0. (25) 

We can combine this new constraint with the original constraint, ^ = 0, and use algebraic techniques to 
obtain an ordinary differential equation for the dynamics of the circuit. 




Figure 5. A parallel RLC circuit, with source: We can use the established rules to write the 
Lagrangian function for the circuit as: 

C = + {jR/2){i R { ~ 1/2) ^j + (L/2) (Vl (0) ) - (1/(2C)) (i c { ~ 1} ) . We also use Kirchhoff 's current law 
to impose a constraint function of VP = ir + i^ + ic ~ is — 0, and we use the principle of the Lagrange 
multiplier to obtain an ordinary differential equation to describe the dynamical behaviour of v(t). 



We can use the Euler-Lagrange terms in Table [5] to write: 

-Ri R ^ 

+W 0) 









V£ = 


5 iL C 






5 ic C _ 





The function of constraint is derived from Kirchhoff 's current law and can be written as ^ 
ic ~ is, and we obtain the gradient of this as 



(26) 



iR + ih + 





" 6 iR * ' 




' 1 " 


V* = 






1 




<5io* 




1 



We can apply the principle of the Lagrange-multiplicr to obtain 



R- 






' 1 " 




YL 




- A- 


1 


= 


i 

c ' 


*c ( " 2) 




1 





for some constant complex number, A G C, which leads to the result 



R ■ i 



,(-!) 



(27) 



(28) 



(29) 



If we apply the constitutive laws for the three devices, and Kirchhoff 's voltage law, three times, then we 
realise that we can interpret A as the common voltage across all three components, A = v. Of course, 
we could have obtained this result using more conventional circuit theory but the point here is that we 
have arrived at differential equations for the system in Figure [5] using purely variational techniques and 
we have been able to model a complete electrical system that includes a dissipative element, R. 
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A ladder filter 

The use of constraints can be a powerful technique, but it does add some extra complication to the 
analysis. It is often possible to make a careful choice of generalised coordinates, and avoid the need for 
constraints. We demonstrate this concept by analysing the ladder circuit in Figure [5] 



If circuits have obvious regularity or symmetry, like the ladder circuit, then it is often possible to choose 
the coordinates in such a way that constraints are automatically obeyed. For the circuit in Figure [51 we 
can choose the state- variables as v\ and V3. The conventional expression for the stored energy in the 
system can be written entirely in terms of these state variables as U = C\ («i) 2 /2 + C2 iv-if /2. It is also 
possible to express all the voltages across the resistors in terms of v\, v 2 and v s , using KirchhofFs voltage 
lav|l. Our decision to use the state-variables as the independent coordinates of the system means that 
the Lagrangian function in Equation [30] takes a simple form and can be written down almost as quickly 
as the circuit can be drawn. Our choice also ensures that the final Eulcr-Lagrangc equations are closely 
related to the state variable model, which could be obtained by using conventional circuit analysis. It 
is also clear, from the symmetry of ladder circuits, that we could extend this Lagrangian technique to 
ladders of arbitrary length and composition, as long as they were composed of components from Table [3] 
The Lagrangian function can be written directly as: 

Rs R2 



vs 







1 — W — I 












+ 






+ 


d = 


— v\ C 3 — 


— v 3 






linear source 


ladder filter 


load 



Figure 6. An electrical ladder-filter circuit: In this circuit, we make a careful choice of 
generalised coordinates, which allows us to avoid the explicit use of functions of constraint. We can use 
the established rules to write down the Lagrangian function for the circuit as shown in Equation [30] We 
apply the rules for the Euler-Lagrange equation to this Lagrangian function to obtain a pair of ordinary 
differential equations that describe the dynamics of this circuit. 



3 (1 \(-i/2) 
. 2 



^3 (0) 



(« 3 ( - 1/2) )' • (30) 



2R L 

This Lagrangian function is completely composed of terms that can be found in Table [3] We can use the 



2 Kirchhoff 's Voltage Law (KVL) does impose constraints on the system. We implicitly use KVL, and apply constraints, 
in the definition of the Lagrangian function, C. This means that we do not need to explicitly use a Lagrange multiplier 
technique. 
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established rules to calculate the variations of the Lagrangian: 

S V1 C = -C lV [ 0) 



Rs 
l 



(31) 



and 



S V3 C 



1 

~R~2 
1 

"rI 



0. 



(32) 



We can differentiate these variations, once with respect to time and rearrange the equations into the 
usual form of a state-variable model of the form 



V = AV + B, 



(33) 



where the time-rate of change is 



and the transition-matrix is 



dt 



Vl 

V3 



R2 J Ci 



+1 

R2C1 



+1 

-R2C3 



o 3 



and the state-vector is 



and the source- vector is 



V = 



Vl 

V3 



B 



RsCi 





Vs- 



There is great advantage in noting the symmetries of any circuit that is being analysed, and choosing 
the generalised coordinates in a consistent manner. For example, we can exploit the symmetry of ladder 
circuits to extend our Lagrangian technique to ladder circuits of arbitrary length. 



A mechatronic problem, the D'Arsonval galvanometer 

One of the great advantages of the Lagrangian approach is that it can be easily used to model devices 
that transduce energy between different forms. For example, an electric motor (or generator) transduces 
energy between electrical energy (in electric and magnetic forms) to, and from, mechanical energy (in 
kinetic and elastic forms). 

The simplest form of electric motor is a piece of wire, moving in a magnetic field. A short element 
of wire, dl, moving in a magnetic field will experience Lorenz force of dF = B ■ i ■ dl, provided that the 
wire and the fields are orthogonal. The modern form of the moving-coil current meter is the result of a 
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long line of development, which includes contributions from many people, including Oersted, Schwciggcr, 
Kelvin, D'Arsonval, Weston, and Ayrton. A typical physical meter is shown in Figure [7J The meter is 
carefully designed to guarantee that the magnetic flux density, B, is orthogonal to the moving wires. We 
can use the Lorenz force on the wire and the radius of the motion of the wire, r, to calculate the rate of 
energy that is transduced per unit angle of motion: 

dU M = (Blr) ■ i ■ dO = $ ■ i ■ d9, (34) 

where <f>o is a parameter that represents the construction of the meter. If we use multiple turns of wire 
then this simply re-scales the parameter, $o, but does not alter the basic model. We can integrate 
Equation IMl to obtain 

U M = $ -i- 0, (35) 

which is the appropriate energy term for the Lagrangian function of a moving, current carrying coil, in a 
magnetic field. We note that the energy that can be transduced is unbounded, if the angle, 9, is allowed 
to increase without bound. Of course, this is normal for a motor. On practice the angle for the meter 
cannot increase outside of the range —90° < 9 < +90° because the meter does not have a commutator. 
Forces would cease and then change direction at the boundaries. 



The moving-coil current meter is shown in Figure [7] The relevant parameters of this physical system 
are, the rotational moment of inertial of the coil (including the needle and the physical supports) J, the 
torsional Hooke's law spring constant, k, the torsional damping constant, \ and the maximum magnetic- 
flux linked by the coil, $o, defined earlier. For a coil in free space the stored energy in the coil is given 
by Um = • i • sin (6^°'), but for the D'Arsonval and Weston style of meter the magnetic field is at right 
angles to a narrow circular air-gap. This makes the magnetic stored energy function more linear with 
respect to 6. We can write Um = $o • i ■ @ i as m Equation 1551 These abstract aspects of the meter and 
their relationships are shown in Figure HI There is only one generalised coordinate for this system, 9 and 
the Lagrangian function for the linearised system can be written as 




+$o^ (0) - (36) 
The last term represents the transduction of energy through the coil. 



We can apply Equation|4]to Equation [36] and obtain the equation of motion for the D'Arsonval meter: 

J9 (2) + x9 {1) + k9 (0) = $ «- (37) 

This example shows that it is possible to model mixed mechanical and electrical (mechatronic) systems 
using Lagrangian techniques. Further, we show that the presence of vicious damping is no obstacle to 
Lagrangian analysis. 



Summary and conclusions 

We have extended the range of applications of Lagrangian analysis, to include non-conservative systems 
that include dissipativc forces. This has been achieved, even though it contradicts many of the accepted 
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Figure 7. Physical layout of the D'Arsonval galvanometer: We model the essential features of 
the D'Arsonval meter as: the rotational moment of inertial of the coil J, the torsional spring constant, 
k, the torsional damping constant, x an( i the maximum magnetic-flux linked by the coil, $o- We use a 
linear model for the stored energy in the coil, U M = $o^ (0) - The La grangian function can be written in 
terms of these fundamental parameters. (Adapted from the Wikimcdia commons.) 



Figure 8. An equivalent electro-mechanical circuit for a D'Arsonval galvanometer: The 

Lagrangian function for this mechatronic system is written in Equation 1371 The current, i, comes from 
an ideal current source, so it is essentially a constraint, rather than an independent coordinate. The last 
term in this Lagrangian function determines the coupling between the electrical and mechanical aspects 
of this system. 

ideas in the current literature. We have also provided a systematic method of applying an extended type 
of Lagrangian analysis to non-conservative mechatronic systems. 

Our approach motivates a number of directions for future work: 

• It is possible to extend Lagrangian techniques to non-linear dissipative systems, such as memristors 
or diodes, using Taylor's theorem, or by using repeated integration by parts, which is how Cauchy 
proved Taylor's theorem. 

• If we could extend fractional calculus of variations to include generalised functions, such as white 
noise, then we could develop a fractional Malliavin calculus. The greater aim is to analyse elec- 
tromechanical systems in the presence of noise. We expect that this would lead to the solution of 





damped 
rol arv 
support 
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the apparent paradoxes of the Penfield motor [B] , and the Davis electromechanical capacitor (37J . 

• Extremal principles can be used to create a number of numerical methods, including the method of 
Rayleigh-Ritz. It is possible to simulate mechatronic systems without having to formulate or solve 
the ordinary differential equations of motion for those systems. It is possible to write down the 
Lagrangian function for a system, specify any constraints, and then find approximate solutions using 
numerical optimisation packages. The authors have had some success with f inconO in Matlab, 
and with sqp() in GNU Octave. Such methods can be iterative, so an approximate solution can 
always be improved, through further iteration 

• It is possible to apply Noether's theorem to determine the "constants of the motion" for quite general 
systems, including systems with dissipative elements and dependent sources. These constants of 
the motion will be generalised forms of momentum and energy. It seems reasonable to suppose that 
the constants of the motion could be used to construct Lyapunov functions, which may then be 
used to make statements about the stability, or otherwise, of mechatronic systems. 

• For noisy electrical systems, with many degrees of freedom, it is of great theoretical interest to 
write down Liouville's theorem, in its most general form. The greater aim here is to understand 
the thermodynamics of electrical systems. 

• It should be possible to create a time-average Lagrangian analysis for switched-mode systems. This 
would be analogous to the time-averaged state-space models of Middlebrook and Cuk [38] . 

In summary, we argue that the generalised Lagrangian functions described in this paper are expected 
to have impact on theoretical and practical applications in electrical and mechanical engineering. 
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